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The ion sphere model introduced long ago by Salpeter is placed in a rigorous theoretical setting. 
The leading corrections to this model for very highly charged but dilute ions in thermal equilibrium 
with a weakly coupled, one-component background plasma are explicitly computed, and the sub- 
leading corrections shown to be negligibly small. Such analytic results for very strong coupling are 
rarely available, and they can serve as benchmarks for testing computer models in this limit. 
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Here we shall describe a plasma configuration that is of 
considerable interest: very dilute "impurity" ions of very 
high charge Z p e, Z p 1, in thermal equilibrium with a 
classical, one-component "background" plasma of charge 
ze and number density n, at temperature T = 1/(3. The 
background plasma is neutralized in the usual way, and 
it is weakly coupled. We use rationalized electrostatic 
units and measure temperature in energy units so that 
the background plasma Debye wave number appears as 
k 2 d = (3 (ze) 2 n. The weak coupling of the background 
plasma is conveyed by g <C 1, where g = (ze) 2 Kr>/(4TrT). 
Although the background plasma coupling to itself is as- 
sumed to be very weak and the impurity ions are assumed 
to be so very dilute that their internal interactions are 
also very small, we shall require that the ionic charge 
Z p is so great that the coupling between the impurity 
ions and the background plasma is large, gZ p 3> 1. This 
strongly coupled system is interesting from a theoreti- 
cal point of view and our results can be used to check 
numerical methods. 

This limit can be solved exactly. The solution is given 
by the ion sphere result presented by Salpeter [l| plus a 
simple smaller correction. This is accomplished by using 
the effective plasma field theory methods advocated by 
Brown and Yaffe |2j. In this field-theory language, the 
old Salpeter result corresponds to the tree approximation 
and our new correction is the one-loop term. 

In usual perturbative expansions, the tree approxima- 
tion provides the first, lowest-order term. Here, on the 
contrary, the tree approximation provides the leading 
term for strong coupling, with the corrections of higher 
order in the inverse coupling. This is the only example 
of which we are aware in which the tree approximation 
yields the strong coupling limit. 

Standard methods express the grand canonical parti- 
tion function in terms of functional integrals. Brown and 
Yaffe 2] do this, introduce an auxiliary scalar electro- 
static potential, and (formally) integrate out the charged 
particle fields to obtain the effective theory. The sad- 
dle point expansion of this form for the grand partition 



function yields a perturbative expansion, with the tree 
approximation providing the lowest-order term. Here, 
on the contrary we express the impurity ion number in 
terms of an effective field theory realized by a functional 
integral. This corresponds to a mixed thermal ensemble 
in which the very dilute impurity ions are represented 
by a canonical ensemble, with the remaining background 
plasma described by a grand canonical ensemble. The 
saddle point of this form of the functional integral in- 
volves a classical field solution driven by a strong point 
charge. 

The result for the impurity ion number reads 

N P =N^ exp |^ (3 . 9 f/3 Z 5/3 

+ ^-y /3 CZ p 2 / 3 + ...-i. 9 Z p | . (1) 

Here ~ exp{/3^ p } is the the number of impurity 

ions defined by the chemical potential \i v in the absence 
of the background plasma; keeping this chemical poten- 
tial fixed, the background plasma alters this number to 
be N p . The final —gZ p /3 term in the exponent is the 
relatively small one-loop correction. The added • • ■ rep- 
resent corrections to the evaluation of the classical action 
that may or may not be significant — if needed, they may 
be obtained numerically. The constant C — 0.8498 

The number correction can be used to construct the 
grand canonical partition function Z for the combined 
system by integrating the generic relation N — ^f- In Z . 
The equation of state is then determined from pV = \nZ. 
To simply bring out the main point, we include here only 
the leading terms, 

pV^^N-Zp^f^Np^T. (2) 

Although the fraction of impurity ions in the plasma 
N p /N may be quite small, there may be a significant 
pressure modification if Z p is very large. 
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The number result Q also directly yields the plasma 
correction to a nuclear fusion rate, since 3] 



r = r, 



n[ 0) n 1+2 

Ni N 2 N (0) 2 



(3) 



where Tc is the nuclear reaction rate for a thermal, 
Maxwell-Boltzmann distribution of the initial (1,2) par- 
ticles in the absence of the background plasma. We use 
the notation 1 + 2 to denote an effective particle that 
carries the charge [Z\ + Z%)e. Thus 



r = r c exp \ — (3. 9 ) 2 / 3 \[Zx+Z 2 
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The first line agrees with the calculation of Salpeter 
Q;the second is new. Again the correction can be large. 

We turn now to sketch the basis for these results. The 
effective field theory expresses 



N (0) 
p Z 



[dx] e 



S[x] 



(5) 



where the effective action is given by 



S[ X ] = |(d 3 r)|^(v X (r; 



Jf3zex(r) 



i(3zex(r) - ipZ p e5(r) x(r) 



(6) 



The normalizing partition function Z is the same func- 
tional integral except that the point source S function 
term is removed from the effective action . The terms 
subtracted from the exponential in the action © remove 
an overall number contribution and account for the ef- 
fect of the rigid neutralizing background. As described 
in Brown and Yaffe one can establish this result by 
expanding the exponential in powers of the number den- 
sity n and performing the resulting Gaussian functional 
integrals to get the usual statistical mechanical form. 

The loop expansion is an expansion about the saddle 
point of the functional integral. At this point, the action 
S[x] is stationary, and the field x obeys the classical field 
equation. The tree approximation is given by evaluating 
S[i<j> c i(r)], where c i(r) obeys the classical field equation 



-V 2 </> cl (r) = zen 



-/3ze0 cJ (r) _ ^ 



ZeS(r) 



(7) 



This equation is of the familiar Debye-Huckle form. We 
have placed it in the context of a systematic perturbative 
expansion in which the error of omitted terms can be 
ascertained. We shall describe the one-loop correction 



automatically produced by our formalism and prove that 
higher-order corrections may be neglected. 

The one-loop correction to this first tree approximation 
is obtained by writing the functional integration variable 
as x( r ) = i<f>c\{ r ) +x'(r) > and expanding the total action 
in Eq. (JSJ to quadratic order in the fluctuating field x' ■ 
Since i<f> c \ obeys the classical field equation, there are no 
linear terms in x' ■ The leading quadratic terms define 
a Gaussian functional integral that produces a Fredholm 
determinant. Hence, to tree plus one-loop order, 
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exp{-S[i<j> c i]} 



(8) 

To solve the classical field equation Q in the large Z p 
limit, we note that <p c \ must vanish asymptotically, hence 
Eq. Q reduces at large distances to the Debye form and 
thus, for |r| large, 



</> c i(r) ~ (const) 



-k\t\ 



(9) 



The coordinate integral of V 2 </> c i vanishes by Gauss' the- 
orem, and from Eq. J7J we obtain the integral constraint 



2/1 



-0e<j> c i(r) 



Z n 



(10) 



For small |r|, the point source driving term in the clas- 
sical field equation dominates, giving the Coulomb po- 
tential solution. Thus we write 



iirr 



(11) 



where £ = nr , and the point driving charge Z p e is now 
conveyed in the boundary condition u(0) — 1 . The large 
r limit © requires that u(£) ~ for large £. 

To compute the action corresponding to the classical 
solution, we must first regularize it and remove the vac- 
uum self energy of the impurity ion. It is not difficult to 
show that this gives, on changing variables to £ = kx, 



[2<M = 
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The classical field equation now appears as 



-9 Z p 



d 2 u{j) 
d£ 2 



= £ 



(12) 



(13) 



according to the variation of Eq. (|12H . 

In our large Z p limit, the short distance form of 
Eq. i|ll[l (multiplied by (3ze) is large over a wide range of 
r|, and exp{— j3ze<j) c \(r)} is quite small there, leading to 
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FIG. 1: Numerical solution for (solid line), the ion sphere 
model uo(£) (short-dashed line), and the first analytic correc- 
tion v\ (long-dashed line). Recall that uo = for £ > £o; here 
£o = 6.694. 



the "ion sphere model" introduced long ago by Salpeter 
. This model makes the step- function approximation 



exp 
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Placing this in the integral constraint (|1(J|) determines 
the ion sphere radius £o = « r o to be given by £q — 3.gZ p . 
Approximating Eq. I|13|) with the replacement Eq. Ijl4(l 
produces a simple differential equation whose solution 
obeying the boundary conditions is 



u (0 



{a/2gz p ) [e - ±c 2 ] 



c<c , 
c>c - 



(15) 



The nature of this "ion-sphere" solution uo(C) together 
with the exact solution obtained by the numerical 
integration of Eq. (| 1 31) , as well as the first correction de- 
scribed below, are displayed in Fig. ^ 

We have appended the subscript to indicate that this 
is the solution for the ion sphere model. Placing this 
solution in the new version 1|12|1 of the action gives 

q y 

[#d] - ^ (3 5 Z P ) 2 / 3 -Z v . (16) 

To find the leading correction, we first write the full 
solution u(£) as u(£) = u (^) + (^o/Z p g) v(£), where i«o(C) 
is the solution l|15|l to the ion sphere model. The exact 
differential equation (|13|) now reads 



d 2 v{j) 



1 

Co 



*(C-Co) 



gZ p u (C) 

e « exp 



Co 

-j «(C) 



(17) 



Since tto(0) = 1 and since the solution must vanish at 
infinity, the proper solution to Eq. (|17|) obeys v(Q) = 0, 
and v(£) — ► for £ — > oo. Some algebra yields 



SVegl^cl] = S a [icj) c i} - 



Co 



I' < 



(18) 



Thus far we have made no approximations. To obtain 
the leading correction to the ion sphere result, we note 
that the factor exp { — {gZ p /^) «o(C)} is vei T small for 
£ < Co, and so it may be evaluated by expanding uo(C) 
about £ = £o- Using Eq. (JT3J, the leading terms yield 



exp 



gZp 
C 



«o(0 
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This approximation is valid for all £ because when £ is 
smaller than £o and our expansion breaks down, the ar- 
gument in the exponent is so large that the exponential 
function essentially vanishes. Since we consider the limit 
in which £o is taken to be very large and the Gaussian 
contribution is very narrow on the scale set by £o, we 
may approximate 



Co) + #(C-Co) . (19) 



Here the delta function accounts for the little piece of area 
that the Gaussian provides near the ion sphere radius. 
One may verify that, with this approximation inserted 
into Eq. fT7f >. the leading solution is given for £ < Co by 
v\ (£) = Ci £ , where c\ is a constant that is yet to be 
determined, while 



£>£o 
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This differential equation is akin to a one-dimensional 
equation of motion of a particle with £ playing the role 
of time, and i>i(C) playing the role of position. Thus there 
is the usual "energy constant of motion" . The integra- 
tion constant is fixed by requiring that vi(C) vanishes at 
infinity. Then choosing the proper root to ensure that 
asymptotically i>i(£) decreases when £ increases gives 



dvi(Q 
d£ 



= -j2 [ e -»i(0 +«!(£)-!] 



(21) 



The different functional forms for Ui(£) in the two re- 
gions C < Co and £ > Co are joined by the continuity 
constraint c\ £o = «i(Co) > and a slope jump to produce 
the 5 function introduced by Eq. Ijl9|) . This requires 



^2 [e-M«o) +Wl (£ )_i] 



,-«i«o) 



MCo) 
Co 



(22) 

For » 1, the second term on the right-hand may be 
neglected, giving «i(Co) = 0.6967 • ■ • . 

We now evaluate the leading correction in the action 
(|18fl . In computing the leading term we can set £ = 
£o in the integral that is linear in «i(C)- The leading 
correction is given by 5 reg [i0ci] — <So[*0ci] + S\ , in which 
5i = -(Co 2 /ff)C, where 
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(23) 
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FIG. 2: Ratios of the modified action S[i(/> c i 
sphere model result lib! [short-dashed line] and the corrected 
ion sphere model (12 H [long-dashed line] to the action 11211 of 
the exact numerical solution as functions of gZ v . 



Here we have omitted a part that is of the negligible 
relative order l/£o- We change variables from £ to v\ 
and use the result (|21[1 for the derivative to get simple 
numerical integrals yielding C = 0.8498 • • ■ . 
In summary, we now find that 



-[Sq + Si]+Z p 



.2/3 



1 + 
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(24) 



The leading correction to the ion sphere model is of rel- 
ative order l/(gZ p ). Fig. [21 displays the exact numerical 
evaluation of the action compared with the ion sphere 
approximation and the corrected ion sphere model. 

The one-loop correction for the background plasma 
with no "impurity" ions present is given by Q 



Det" 1 / 2 [- V 2 + k 2 ] = cxp U (d 3 r) ^- J 



(25) 



We assume that the charge Z p of the "impurity" ions is 
so large that not only Z p ^> 1, but also gZ p ^> 1 as well, 
even though we require that g <C 1. Then ktq 3> 1, and 
the ion sphere radius tq is large in comparison to the 
characteristic distance scale for spatial variation in the 
background plasma, the Debye length In this case, 
the term k 2 exp {— f3ze(j)(r)} in the one- loop determinant 
Det~ 1/2 [-V 2 + k 2 e -^ e *=i] i n Eq. (JgJ can be treated as 
being very slowly varying — essentially a constant - 
except when it appears in a final volume integral akin to 
that in Eq. Therefore, for very strong coupling, 
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= exp j- ^- J (d 3 r) [1 - cxp {-pz 
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exp 
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(26) 



where in the second equality we have used the ion sphere 
model that gives the leading term for large Z p . 



This result is physically obvious. The impurity ion 
carves out a hole of radius rrj in the original, background 
plasma. The original, background plasma is unchanged 
outside this hole. The corrections that smooth out the 
sharp boundaries in this picture only produce higher- 
order terms. The original, background plasma had a 
vanishing electrostatic potential everywhere, and the po- 
tential in the ion sphere picture vanishes outside the 
sphere of radius rrj. The grand potential of the back- 
ground plasma is now reduced by the amount originally 
contained within the sphere of radius rrj, and this is ex- 
actly what is stated to one-loop order in Ea. (|2"6*|) . 

This argument carries on to the higher loop terms as 
well. As shown in detail in the paper of Brown and Yaffe 
Q, ro-loop terms in the expansion of the background 
plasma partition function with no impurities present in- 
volve a factor of k 2 n n which combines with other factors 
to give dimensionless terms of the form <?™ _1 J (d 3 r) k 3 . 
With the very high Z p impurity ions present, each factor 
of k is accompanied by exp{ — (l/2)/3ze cf> c \(r)} whose spa- 
tial variation can be neglected except in the final volume 
integral. Thus, an n-loop term is of order 



(d 3 r) 



1 — exp < — 



„n-l „3„3 



■Pze<f>d(r) 



(27) 



Since we assume that g is sufficiently small that g 2 Z p <C 
1 (even though gZ p 3> 1), all the higher loop terms may 
be neglected |5J. 

We have now established the results quoted above. 

We thank Hugh E. DeWitt and Lawrence G. Yaffe for 
providing constructive comments on preliminary versions 
of this work. 
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